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Abstract 



We investigate the low-energy spectrum of excitations of a compressible elec- 
tron liquid in a strong magnetic field. These excitations are localized at the 
periphery of the system. The analysis of a realistic model of a smooth edge 
yields new branches of acoustic excitation spectrum in addition to the well 
known edge magnetoplasmon mode. The velocities are found and the observ- 
ability conditions are established for the new modes. 
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The dispersion relation for plasmons in a non-restricted two-dimensional electron liquid 
is well known to have a form lo oc k 1 ^ 2 |l]-f|. If the liquid has a boundary, an edge mode 
appears in addition to these bulk excitations. The spectra of the edge and bulk modes differ 
from each other only by a numerical factor [|J. A magnetic field applied perpendicularly 
to the plane of the liquid changes the plasmon spectrum drastically. The spectrum of the 
bulk mode acquires a gap of the width equal to the cyclotron frequency uj c . The only known 
gapless mode existing in the presence of the magnetic field propagates along the boundary 
0-[7|]. The "chirality" of this edge magnetoplasmon determined by the direction of the 
magnetic field (i. e., by the sign of the Hall conductivity a xy ) was demonstrated explicitly 
in the time-domain experiments ||. 

The solved theoretical models of the edge modes assumed a sharp electron density profile 
at the boundary 0-0], i. e. width of the boundary strip was assumed to be infinitesimal. 
The existence of only a single branch of the edge magnetoplasmons follows directly from this 
assumption. For a realistic shape of a potential confining the electron liquid, the density 
profile is smooth at the boundary [fl-fOll- The results of Refs. f§-[7| can be extended on 
this case only under the assumption that the current and charge oscillations forming the 
magnetoplasmon wave are homogeneous across the boundary strip. However, the latter 
condition is excessively restrictive. We demonstrate in this paper the existence of other 
sound-like modes propagating along the edge. The current for each of these modes alternates 
across the boundary strip, and therefore the new branches could not be predicted on the 
basis of a "sharp" boundary model. 

Below we present an exactly solvable model correctly describing all the edge excitations 
in the strong magnetic field limit. We obtain also the values of the oscillator strengths and 
the damping of these modes. The new branches become robust and are not destroyed by a 
finite relaxation time in the achievable region of relatively short wavelengths. 

The dynamics of the compressible electron liquid is governed by Euler equation and the 
continuity equation linearized in the velocity of the liquid v (p, t) and in the deviation of 
the concentration 5n (p,t) from its equilibrium value no (p): 

v + co c (zxv)- £v p /rf 2 Pl ^ = 0. (!) 

5h + V p (n v) = 0. (2) 

Here p is radius-vector in the plane XY of the two-dimensional electron liquid, z is the 
unit vector along Z-axis, and e is the dielectric constant. The last term in ([!]) represents 
Coulomb interaction [I2|. 



In the following we assume that the electron liquid is homogeneous in y direction and 
occupies half-plane x > 0. Since the system is translationally invariant in the y direction, 
we will seek the solution of Eqs. ([!]), (0) in the form 

v = exp (iky — iuot) w(x); 
Sn(x, y) = exp (iky — iut) f(x). (3) 

Substituting Eqs. (§) into the system (p]), (0) and eliminating w(x), we find an integral 
equation for f(x): 
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d 2 d k 

k 2 n - n o-j^ ~ n' — + -u c n' \ / K (\k\\x - x t \) f(xi)dxi = 0, 



u 



(4) 



where n' Q = dn /dx and i^o(^) is the modified Bessel function. Homogeneous equation 
(B) comprises the eigenvalue problem that determines the spectrum of edge excitations 
u)j(k). The spectrum is controlled by the parameters of the problem: by the magnetic field 
determining u c , and by the concentration profile n (x). The latter depends on a particular 
type of the confining potential [§,[n|] • We are interested in the low-frequency, long wavelength 
modes, and this allows us to neglect the terms proportional to u 2 and k 2 in Eq. (£|). 

Futher simplifications are possible in the case of a strong magnetic field. Keeping only 
the terms proportional to u 2 and u c , and introducing a new function 

instead of f(x), we find from @ the following reduced equation: 



g(x) = X [ K (\k\\x - x 1 \)-J ( ^- ( ^-g(x 1 )dx 1 , (6) 
Jo n V dx ax\ 
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u = — - 



X emu, 



k. (7) 



Here n = uq(x — > oo) is the density of the homogeneous electron liquid far from the bound- 
ary. One can estimate from Eqs. (j|), (0) the typical value of \u/k\ to be of the order of 
fie 2 /mu c . It follows also from (^j), (|]) that the charge distribution f(x) in the wave is local- 
ized mainly within the region where dn /dx is large, i. e. within the boundary strip of width 
a. Now we can establish the validity criterion of the strong magnetic field approximation. 
The neglected in Eq. (f|) terms oc n/a are smaller than the main terms oc mu 2 /e 2 if the 
condition 

» — (8) 

ema 

is satisfied. For realistic parameters of the two dimensional electron system formed in a 
GaAs heterostructure |TD|, n ~ and a ~ 10a b, the latter condition is equivalent to a 



rather weak restriction on the filling factor, v < 10. 

Eq. (^j) is the integral equation of Fredholm type. Its kernel is symmetric and positively 
defined, hence, all the eigenvalues A are real and positive. If one makes an approximation 
-Kod^H^ - x\\) fa In (l/|fca|) leading to the degeneracy of the kernel, then only a single finite 
eigenvalue A exists. This eigenvalue corresponds to the known magnetoplasmon mode [|J. 
The actual kernel in is non-degenerate, however, and thus there are many edge modes. 

The eigenvalue problem (^]) can not be solved analytically for an arbitrary distribution 
uq(x). Below we present a model for the density profile : 

2 fx~ 

n (x) = — narctanW— , x > 0, (9) 
7r V a 
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that allows a complete analytical solution of the problem. This model describes correctly 
the asymptotic behavior of the density formed by an electrostatic confinement, reproducing 
the characteristic |9|-|TT| \fx- singularity at x — > 0. 



Proposed model allows us to solve the eigenvalue problem @ using an expansion of g(x) 
in a form 



g(x) = (W+ay^l^ (v^j ' (10) 

where T n (£) = cos(narccos£) are the Chebyshev polynomials |^3| . Substitution of (|10|) into 
(§) leads to the following system of equations for coefficients gj of the expansion: 




(11) 

When deriving Eqs. flTT|), we used the approximation Ko(kx) ~ In {2e~ 1 /\kx\) which is valid 
in the long-wavelength limit, \ka\ <C 1; here 7 ~ 0.577... is the Euler constant. The system 
(PI) leads directly to the following transcendent equation for the eigenvalues: 

x + 2 *< 1 - A) = ln (lS)- (12) 



Here ^(1 — A) is the digamma function [Tj| that has simple poles at A = 1, 2, . . .. Because 



|fca| <C 1, the solutions of Eq. fll2"|) are close to the points A = 0, 1,2, . . . where the left- 
hand side of this equation is singular. The smallest root of Eq. (|l^) belongs to the region 
A< 1. Expanding the left-hand side of this equation in power series in A and retaining only 
the two leading terms of the expansion, we find the spectrum of the conventional Q edge 
magnetoplasmon mode: 



u (k) = -2 In ^— ——k. (13) 




Other roots A > 1 are close to the poles of the digamma function in Eq. (Hp. It is these 



roots that determine the new branches of the edge excitations with acoustic spectrum: 

2ne 2 

LVj{k) = -sjk, sj = r , j = 1,2,.... (14) 

The difference between the acoustic modes and the "usual" plasmon (J = 0) originates 
in the structure of charge distributions associated with these waves. In the usual plasmon 
wave, charge does not oscillate across the boundary strip, whereas in the acoustic mode j 
charge oscillates j + 1 times in x direction (see Fig. 1), the average density being smaller 
by a factor of \j ln(|fca|) | _1 than the oscillations amplitude for each mode with j 7^ 0. The 
potential energies Uj produced by the charge distribution types depicted in Fig. 1 can be 
easily estimated. For the same characteristic amplitudes of charge density perturbation in 
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all the waves, we find the ratio Uo/Uj ~ j\ ln(|A;a|)|. This explains the difference between 
the spectra (p~3|) and ([14]), as energies Uj provide the restoring forces for the modes. For 
the higher harmonics j ' ^> 1 the latter considerations are obviously independent on the 
particular density profile that allows one to expect certain universality of the spectrum 
(0). Indeed, for an arbitrary profile no(x), the function 

/ dn \ 1/2 ( .n {x)\ 

may be used at j 1 as an asymptotic solution of the eigenvalue problem (§). In the case 
of the profile derived in Ref. ||, we find the corrections to the velocities Sj (see Eq. (|T4|)) to 
be of the order of 1/j 2 . 

The observability of the acoustic modes requires sufficiently large oscillator strengths 
and rather slow decay for these modes. To begin with, we evaluate the oscillator strengths 
S"^(k) for all the modes; here a, (3 = x, y denote the polarization of the applied AC electric 
field, E a (y, t) = E a exp(iky — iut), and j = 0, 1, 2, . . . is the mode number. The power P 
absorbed from the AC field within the unit length of the boundary is related to the oscillator 
strengths of the different modes by: 



P = ^-Re f n vE*d 2 p 



-I oo 

= -T,Sf(k)E a E;5(u J -uj J (k)). (16) 
z j=o 

Here L is the length of the boundary, and velocity v is the linear response to the external 
electric field E. To calculate v , one has to add the term eE/m to the right-hand side of the 
equation of motion ([[]). The use of Eqs. (P, (H) - (0) and (|16D allows to express Sj(k) in 



terms of the eigenfunctions gj(x). We present here the results for the diagonal components 
of the oscillator strengths: 

= ^M^ Fm (17) 




functions gj(x) here being normalized by condition 

J dxg%x) = 1. (18) 

It is obvious from Eq. (|J) that Fj is proportional to the square of the average charge mode j 
bears. As was mentioned already, this charge in the acoustic modes is parametrically smaller 
than the one in the usual edge magnetoplasmon mode. Therefore, the interaction of AC field 
polarized along the boundary with the acoustic modes is much weaker than the interaction 
with the magnetoplasmon. An explicit calculation in the framework of our model gives: 
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= 1 ( ne 2 \\ / 27r[ln(lfeal)l , j = 0; 
j eVmwJ \ 47r|ln(|fca|)r 2 j~ 3 , j > 1. ^ 

AC electric field applied perpendicular to the boundary interacts with the x component of 
the dipolar moment of the modes. These moments, and correspondingly factors Fj are of 
the same order of magnitude for all the modes. Therefore, the difference in absorption for 
this polarization is only due to the difference in the mode frequencies: 




27r 3 |ln(|H)r\ j = 0] 

47r 3 |m(M)r 2 r 1 , j>i- 



aiwi^m-a ,-! , ^ •. ( 20 ) 



It is interesting to notice that the absorption anisotropies SJ X / Sj V — 1 are of the opposite 
signs for the usual plasmon and for the acoustic modes respectively. 

To estimate the decay rates of different edge modes, we include a phenomenological 
relaxation time r into the equation of motion ([!]) by the substitution v —>■ v + v/r. After 
such a modification the energy of the mode 6j becomes time dependent, 6j(t) oc exp(— 2t/r,). 
The latter relation allows one to define the relaxation rates 1/tj. For small dissipation the 
result can be expressed in terms of the unperturbed eigenmodes gj(x) : 

^ k u ixna(x) (ijA\ (21) 



tj uj c rk J \dx 



gj(x) being normalized by the condition (0). As one can easily see from Eq. (pl|) , the 
relaxation rate increases with the mode number because of oscillatory behavior of the eigen- 
functions. We find for the plasmon and acoustic modes: 

^ = H-orx{»' 2 ^: = , ; ; (22) 

where /3i = 6/5, /?2 = 60/53,. . . , = 1.128 are slowly varying with j numerical factors. 

Solutions of Eq. (||) and perturbative results ( |22"D obtained above are applicable as long 
as dissipation is small enough so that it does not affect the charge distribution in the eigen- 
modes. The characteristic length of the charge spreading Z w caused by dissipation is inversely 
proportional to frequency, = e 2 n/(emuJU!^T). The smallness of the redistribution requires 
1^ to be shorter than the characteristic length scale ~ a/(J + 1) for spatial variations of the 
eigenfunctions gj(x). The latter condition imposes different restrictions for the plasmon and 
acoustic modes; with the help of (|13"1), QUI) , we find: \ka\ ln(l/|/ca|) ^ \/{uj c t) for j = 
and \ka\ > j 2 /(uj c t) for j ^ 1. At smaller wavevectors, the results of Volkov and Mikhailov 
Q for the spectrum and decay rate are applicable [Q, whereas the acoustic modes are 



over-damped. The region of observability of the new branches is shown on Fig. 2. 

The observability condition for the new modes is quite restrictive. Indeed, as it follows 
from (|22|), the product uj c t must be at least greater than l/|/ca| for the first acoustic mode to 
be observed. If the characteristic values of k are determined by a sample perimeter (which 
is typically ~ 1 cm |T5|), the condition (^) can not be satisfied even for the mobility of 
10 6 cm 2 /V sec. The technique using a metallic grating coupler |16[ appears to be more 



promissing. In such an experiment the wave vector k = 2n /d is determined by the grating 
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period d that may be made |16| of the order of 1 fim. For the typical width of the boundary 
strip ~ 2000A this implies the condition uo c r > 1. 

In conclusion, we found the new low-frequency excitations propagating along the edge 
of a two-dimensional electron liquid in the presence of a magnetic field. These new modes 
have acoustic spectra with velocities inversely proportional to the mode indices. At a given 
wavevector k, the frequencies of acoustic modes are lower than that of a conventional edge 
magnetoplasmon by a factor 1 / 1 In fca | , where a is the width of the boundary strip. The 
oscillator strengths and decay times for the acoustic branches with low indices differ from 
the corresponding values for the magnetoplasmon by powers of the same small parameter. 
The logarithmic function appears due to the long-range nature of Coulomb interaction. If 
the electron system is confined by gate-induced potentials, the differences in the mentioned 
parameters for the acoustic modes and the conventional magnetoplasmon become less signif- 
icant: the logarithmic function should be replaced by a factor of the order of unity because 
of the screening of Coulomb interaction by the metallic gates. 

The authors are grateful to H. U. Baranger, K. A. Matveev, V. I. Perel and B. I. 
Shklovskii for helpful discussions and comments. This work was supported by NSF Grants 
DMR-9117341 and DMR-9020587. 
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FIGURES 



FIG. 1. Characteristic charge distributions for: (a) the edge magnetoplasmon mode, j = 0; 
and (b) the edge acoustic mode with j = 2. Charge patterns shown in the figure move along the 
Y-axis according Eq. (||). 

FIG. 2. The first three branches of the edge excitations; = 2ne 2 / '(emco c a) . The dashed line 
separates the regions of strong damping (below the line) and weak damping (above the line). In 
the latter region the acoustic edge modes become observable. 
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